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ABSTRACT 



Context. Clusters of galaxies are effective gravitational lenses able to magnify background galaxies and making it possible to probe 
the fainter part of the galaxy population. Submillimeter galaxies, which are believed to be star-forming galaxies at typical redshifts of 
2 to 3, are a major contaminant to the extended Sunyaev-Zeldovich (SZ) signal of galaxy clusters. For a proper quantification of the 
SZ signal the contribution of submillimeter galaxies needs to be quantified. 

Aims. The aims of this study are to identify submillimeter sources in the field of the Bullet Cluster (IE 0657-56), a massive cluster 
of galaxies at z ^ 0.3, measure their flux densities at 870 yum, and search for counterparts at other wavelengths to constrain their 
properties. 

Methods. We carried out deep observations of the submillimeter continuum emission at 870 /jm using the Large APEX BOlometer 
CAmera (LABOCA) on the Atacama Pathfinder Experiment (APEX) telescope. Several numerical techniques were used to quantify 
the noise properties of the data and extract sources. 

Results. In total, seventeen sources were found. Thirteen of them lie in the central 10 arcminutes of the map, which has a pixel 
sensitivity of 1.2 mJy per 22" beam. After correction for flux boosting and gravitational tensing, the number counts are consistent 
with published submm measurements. Nine of the sources have infrared counterparts in Spitzer maps. The strongest submm detection 
coincides with a source previously reported at other wavelengths, at an estimated redshift z - 2.7. If the submm flux arises from two 
images of a galaxy magnified by a total factor of 75, as models have suggested, its intrinsic flux would be around 0.6 mJy, consistent 
with an intrinsic luminosity below 10'~Lo. 



Key words. Galaxies: individual: MMJ065837-5557.0 
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1. Introduction 

The large concentrations of mass (up to lO'^M©) on angular 
scales of a few arcminutes in galaxy clusters act as natural grav- 
itational lenses capable of magnifying background galaxies that 
would be too dim to be detectable otherwise. In the mm and 
. ^ submm wavebands, lensing by galaxy clusters makes it possible 
K> I to probe the fainter part of the brightness distribution of the so- 
d called submillimeter galaxies (SMGs), w hich are beli eved to be 
d ^ dusty high-redshift star-forming galaxies (lBlainlll997l) . 

Pioneering observations of SMGs at 450 and 850 /im were 
done using SCUBA on the James Clerk Maxwell Telescope. The 
first observations toward two massive clusters at z ^ 0.35 re- 
sulted in the detection of a total o f six sources above the noise 
level of 2 mJy/beam at 850 yum dSmail et al.lll997h . The au- 
thors estimated the surface density of the sources to be three 
orders of magnitude larger than the expectation from a non- 
evolving model using the local IRAS 60 /^m luminosity func- 
tion. Those observations provided evidence for the presence of 
a large number of actively star-forming galaxies at high red- 
shift, which might be the counterparts of the luminous and ultra- 
luminous infrared galaxi es observed in the local universe (e.g. 
ISanders & Mirabeff996l) . 



At redshifts beyond one, the flux density of a redshifted 
infrared-luminous galaxy is largely redshift-independent : its de- 
crease with an increasing distance is compensated by the steep 
rise in the mm and submm dtie to t he redshifted spectral energy 
distribution ( lBlain&Longaijl993h . During the last decade, sev- 
eral hundreds of SMGs have b een discovered using bolome- 
ter arrays, mostly S CUBA (e.g. lBlainlll998l: iBorvs et al.ll200l 
ICoppin et al] |2006|) . and more recently MA MBO at 12 mm 
and AzTEC at 1.1 mm jBertoldi et al. 2007; Scott etal.ll200l 
lAustermann et al.ll2010l) . The mm/submm galaxy population is 
th e subject of many multi-wavelength studies (see the review 
bv lBlain et al]|2002l) . Th e median redshift o f SMGs with known 
redshifts is around 2-3 (ISmail et al.ll2000h . 

So far, only a handful of massive galaxy clusters have been 
mapped in the submm, and most of them are clusters in the 
northern hemisphere observed with SCUBA. Observation of 
seven massive clusters with a sensitivity of 2 mJy/beam pro- 
vided a catalogue o f 17 submm sourc es brighter than the 50% 
completeness limit (ISmail et al.lll998h . Nine other cluster fields 
in the redshift range 0.2 - 0.8 were observed with a similar 
sensitivity, resulting in the detection of 17 new submm sources 
(Chapman et al. 2002). Deeper observations with a 3-sigma limit 
of 1.5-2 mJy/beam of three massive clusters probed the sub- 
mJy number counts, because of the gravitational magnification 
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Fig. 1. Left: layout of the bolometers on the LABOCA array. Each circle represents one bolometer that was active at the time of the 
observations. The two rectangles illustrate the sizes of the two raster+spiral scanning pattern for one bolometer (solid line-small 
pattern, dashed line-large pattern). Right: movement pattern for one bolometer during the large scanning pattern. This pattern has a 
2x2 raster setup with 100" between the raster points (the patterns are described in detail in Sect. 12. lb . The smaller scanning pattern 
is similar, but because the distance between raster points is considerably smaller, it is harder to illustrate. A comparison between the 
two patterns can instead be made from the two rectangles in the left panel. For clarity, we plot also the dashed rounded rectangle of 
the left panel in the righ panel. 



of the clusters. (ICowie et al.ll2002h . iKnudsen et"aD (l2008h tar- 
geted twelve clusters and the New Technology Telescope Deep 
Field. They detected 59 sources (some of them being multiple 
images of the same galaxy), and determined that seven of them 
have sub-mJy lensing-corrected flux densities. 

The LABOCA bolometer camera on APEX has been used 
to survey the 870 yum emission in a protocluster at z ^ 2.4 
(iBeelen et al.l2008h and the Ext e nded C handra Deep Field South 
IWeiB et al.1 ( l2009l) . iNord et al.l (|2009|) observed the Sunyaev- 
Zeldovich (SZ) increment toward the massive cluster Abell 2163 
and noted one bright point source. This is a good candidate for 
an SMG lensed by the cluster 

The Bullet Cluster (IE 0657-56) at z ^ 0.3 is one of the most 
massive galaxy clusters know n to date (see iMarkevitch et alj 
120021: ISoringel & Farrai]l2()07l) . A bright millimeter source was 
recently discovered in the Bullet Cluster field and identified 
as the len sed image of a ba ckground galaxy at a redshift of 
about 2.7 dWilson et alJl2008l hereafter W08). The source hap- 
pens to lie close to a critical line of the lens, causing a large 
flux amplification. The observations were performed with the 
AzTEC bolometer camera on the 10-meter ASTE telescope in 
the Atacama desert in Chile, which provides an angular resolu- 
tion of 30" at 1 . 1 mm wavelength. A doubly lensed source at the 
same location had been previously identified in Hubble Space 
Telescope {HS T) maps and was used, together with other multi- 
ply lensed galaxies and a large number of weakly lensed sources, 
to obtain a calibra ted map of the proj ected mas s distribution of 
theBullet Cluster ( iBradag et alj|2006l) . Recentlv. lOonzalez et al] 
(I2009h identified a third image by analyzing maps in the op- 
tical (HST), and in the near- and mid-infrared (Magellan and 
Spitzer). By fitting the Spectral Energy Distribution (SED) of a 
starburst galaxy to the observations, they also inferred a redshift 
of about 2.7. Their lensing model gave a magnification of 10-50 
for the three images. 

In this paper we present results of observations of the 
Bullet Cluster field at a wavelength of 870 /im using the 



LABOCA bolometer camera. At that wavelength, the emis- 
sion is a combination of extended signal due to the Sunyaev- 
Zeldovich effect by the hot intracluster gas and of point 
sources, which are potential high-redshift star-forming galaxies. 
Recovery of the extended SZ signal from the LABOCA data (the 
SZ increment) requires a different data reduction and will be pre- 
sented in a subsequent paper The SZ decrement from the Bullet 
Cluster has been mapped by e.g. iHalverson et al. (2009) with 
the APEX-SZ instrument, operating at 2 mm. This paper is or- 
ganized as follows: the observations are presented in Sect. 12.11 
and the data reduction in Sect. [3l the results are presented and 
discussed in Sects. 4 and|5] 

Throughout the paper, we adopt the following cosmological 
parameters: a Hubble constant Hq = 70 km s"' Mpc"', a matter 
density parameter Qq = 0.3, and a dark energy density parameter 
Qao = 0.7. The redshift z = 0.296 of the Bullet Cluster corre- 
sponds to an angular-diameter distance of 910 Mpc and a scale 
of 4.41 kpc/arcsec. 

2. Observations 

2.1. Submillimeter 

The observation^ were carried out in September, October and 
November 2007 using LABO CA (Large APEX BOlometer 
CAmera. ISiringo et al.l |2009|) on the APEX telescop^ 
(iGiisten et al.l l2006h . LABOCA is a 295 -element receiver 
operating at a central frequency of 345 GHz with a bandwidth 
of 60 GHz. At the time of our observations, about 250 bolome- 
ters were used. The mean point-source sensitivity of those 



' Swedish program ID O-079.F-9304A 
and ESO program ID E-380.A-3036A 

- This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and Onsala Space Observatory. 
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bolometers was 78 mJy s ' . The angular resolution was 19.5" 
and the field-of-view was 11.4'. The layout of the bolometer 
array is illustrated in Fig. \T\ A total of 25 hours of observing 
time was spent, including pointing and calibration. The weather 
conditions were varying, with an amount of precipitable water 
vapor ranging from 0.5 to 2.0 mm. 

LABOCA uses feed horn antennas and their physical size 
limits the spacing of the bolometers on the array. LABOCA is 
therefore not fully sampling the sky, so the telescope has to be 
moved to sample the sky such that the resulting map meets the 
Nyquist criterion. The scanning pattern does not only "fill the 
gaps" between bolometers, it also modulates the astronomical 
signal into a range of spatial frequencies which facilitates filter- 
ing of 1 / /-type noise (instrumental and sky noise). 



2.1.1. Scanning patterns 

We used two different scanning patterns, as illustrated in Fig.[T] 
Both are a combination of Archimedian spirals with a duration 
of 35 seconds, centered on a four-point raster 

- During the first observing session, we used a compact scan- 
ning pattern: the four points were separated by 27" in azimuth 
and elevation, each point marking the center of a spiral with a 
minimum radius of Ro - 18" and winding out with a radial 
speed of ^ = 2.25"i"' and an angular speed of ^ - 90deg 
The spirals thus ranged from 18" to about 97" in radius, with a 
scanning speed between 0.5'i"' and 2.5'^"'. 

- During the second observing session, a larger scanning pat- 
tern was used to facilitate the retrieval of the extended SZ signal. 
The four raster points were separated by 100" in azimuth and el- 
evation, and the spirals had a minimum radius /?() - 120" and 
radial and angular speeds of ^ = 1.25"i"', = 90degi"'. The 
spirals thus ranged from 120 to 164" in radius, with a scanning 
speed between 3 and 4'i"'. 



2.1.2. Pointing, focus and calibration 

The pointing accuracy was checked by repeated observations 
of the nearby source PKS 0537-441. This is a variable source; 
during our observations, its mean flux density was ~ 3 Jy. The 
source was scanned in a tight spiral and the data were reduced 
and made into a map using the BoA software (see Sect. [3]l. A 
two-dimensional Gaussian was fitted to the pointing source, and 
the telescope's pointing was updated using offsets from the fit. 
The pointing was stable within 3". 

The focus was checked at least twice during every observ- 
ing session by observing a planet (Venus, Saturn or Mars). The 
subreflector was moved in small increments in each of the three 
cartesian directions while the telescope tracked the source. The 
optimum focus position in each direction was determined by fit- 
ting a curve to the observed points, and the subreflector was fi- 
nally moved to the position corresponding to the maximum of 
the curve. 

The absolute flux calibration of L ABOCA is supposed to be 
accurate to 10% (ISiringo et al.ll2009l) . We verified this by daily 
observations of Uranus. 

The atmospheric attenuation was determined from continu- 
ous scans in elevation at a fixed azimuth ("sk ydips"), and from 
radiometer measurements dSiringo et al.ll2009l) . 



2.2. Infrared 

In Sect. 15. II we describe a search for infrared counterparts to the 
detected submm sources in Spitzer maps. We now describe the 
data that was used for that comparison. 

Spitzer IRAC and MIPS of the Bullet Cluster field were ac- 
quired from the Spitzer data archive. Both the IRAC and MIPS 
data were taken under program ID 40593 (PI: Gonzalez). The 
Spitzer data overlap with most of the the region observed with 
LABOCA, and out of the detected submm sources only one 
source lacks Spitzer coverage. 

The IRAC and MIPS maps were processed by, respectively, 
version 18.7 and 18.12 of the SSC pipeline. We started by vi- 
sually inspecting the resulting pbcd (post basic calibrated data) 
mosaics, and found that the IRAC maps had no apparent arte- 
facts, but that the MIPS map had clear signs of "dark latents", as 
described in the MIPS data handbook. We therefore reprocessed 
the basic calibrated data (bed) using MOPEX version 18.3.3 the 
script flatfield.pl to self-calibrate the bed data. The bed's 
were then mosaiced using mopex.pl. The reprocessed MIPS 
map shows no signs of artefacts due to "latents". 

Properties of the acquired Spitzer data are summarized in 
table \T\ We list there the median integration time per pixel per 
map, the sensitivity and the angular size of each map. The sensi- 
tivity is estimated from the final combined mosaics, masking out 
all pixels brighter than 10 times the median pixel value in each 
map, and then calculating the standard deviation of the remain- 
ing pixels. The listed sensitivity values are 3cr, and are similar 
to those obtained from the Spitzer Science Center "Sensitivity - 
Performance Estimation Tool'l^. 

Table 1. Properties of the Spitzer data used in this study. 
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3. Data reduction 

In this section, we describe the steps followed to produce a fully 
calibrated map from the raw data, which come in the form of 
time-streams containing the voltage read-outs of each bolometer 
as a function of time. We have used two data reduction soft- 
wares: MinicrustQ, written originally for the SHARC bolome- 
ter array and adapted to handle LABOCA data (Kovacs 2008j), 
and Z?o4Il developed in Bonn (Schuller et al., in prep). In gen- 
eral, the maps produced by both pipelines were in good agree- 
ment in terms of number and characteristics of sources; however, 
since the Mini crush map showed a lower level of large-scale 
noise, we used that software for the analysis presented in this 
paper. 



^ http : // SSC . spitzer . caltech . edu/ tools/ senspet/ 
Mini crush can be downloaded from 
http : //www. submm. caltech. edu/~sharc/crush/ download. htm 

' BoA can be downloaded from 
http : //www. apex-telescope .org/bolometer/laboca/boa/ 
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The data consist of a total of 185 eight-minute-long scans on 
the cluster, plus pointing and calibration observations. Each scan 
is contained in a separate MBFITS-file. 

First, we flagged blind bolometers to exclude them from the 
rest of the data analysis. Then, we corrected for differences in 
sensitivity of individual bolometers. That information was ex- 
tracted from "beam-maps", which are fully sampled maps where 
each bolometer has scanned a bright planet. 

For each scan, the zenith opacity calculated from the ra- 
diometer and sky dip measurements was used to calibrate the data 
for elevation-dependent opacity variations. We also flagged data 
taken during periods of high telescope speeds and accelerations. 

3.1. Removal of correlated noise in the time-streams 

The most critical and challenging task of the data analysis is 
to extract the true astronomical signal from the measurements, 
which are contaminated by noise from various sources. One 
component of the noise arises in the electronic systems, such 
as the readouts of the bolometer array; its spectrum is of the 1 // 
type where / is the frequency. Even more important is the con- 
tribution of the Earth's atmosphere, which also has a l//-type 
spectrum, but shows both spatial and temporal fluctuations. In 
the 870 yL/m atmospheric window in which LABOCA operates, 
typical zenith opacities at the APEX site are of the order of 0. 1 
to 0.2. The atmosphere is thus largely transparent to cosmic sig- 
nals; but the amplitude of the atmospheric signal can be as high 
as 10^ times that of the astronomical signal of interest. Because 
the atmospheric signal (or noise) is correlated across the bolome- 
ter array, it can be estimated and removed from the time-streams. 
For point-source observations, correlated sky noise can be partly 
removed by filtering low spatial frequencies (or large angular 
scales), where the l/f noise is most severe. 

The software Mini crush takes the following approach to re- 
move correlated noise from the time-streams: a correlated noise 
component is modeled as a common signal in the time-streams 
of the different bolometers, scaled by a gain factor which de- 
pends on each bolometer A function is minimized by tak- 
ing its derivative with respect to the modeled signal. The time- 
streams of all the bolometers are considered for a certain scan. 
When the fit has been performed, the estimated correlated sig- 
nal is removed from the data. The process of estimation and re- 
moval of the correlated signal is carried out several times. The 
uncertainties on the modeled correlated signal are estimated by 
calculating the changes in the estimated signal when the x~ has 
increased by 1 from its minimum value. The ideas implemented 
in Minicrush have been described in detail in Kovacs (2008). 

Since we were interested in compact sources, we used the 
option -deep in Mini crush. 

3.2. Map-making 

3.2.1. From time-streams to maps 

The data consist of a collection of time-streams from each func- 
tioning bolometer, cleaned from atmospheric and instrumental 
noise. The data in each time-stream of a scan are used to create a 
map of that scan, taking into account the scanning pattern of the 
telescope and using a nearest-pixel mapping algorithm. We refer 
to those maps as scan-maps. For each scan-map, a corresponding 
noise map is generated by adding in quadrature the noise levels 
of the bolometers that hit a certain pixel in the map. The noise 
maps are used to weigh the individual scan maps when they are 
co-added, which is the final step of the map-making process. The 
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Fig. 2. Color image of the 870 fiva signal-to-noise map. The 
map has been filtered to remove extended signal, including that 
due to the Sunyaev-Zeldovich effect. Several sources are visible. 
The brightest coincides with the millimeter source discovered by 
W08 and identified as a galaxy at redshift z - 2.1 and strongly 
lensed by the Bullet Cluster The contours refer to the noise level 
at 1.3, 2.1, 2.7, 4.1 and 6.8 mJy/beam. The signal-to-noise rep- 
resentation supresses the high noise levels in the outer parts of 
the map. They are instead indicated in the contours of the noise 
map. 



individual noise maps are also coadded into a final noise map. 
This noise map reflects the pixel rms in the data. As the noise 
map only accounts for the relative weights between map pixel, 
we then reseated the noise map. In Sect. 14.11 we show that the 
signal-to-noise map has a Gaussian distribution of pixel values. 
We scaled the noise map in order for the signal-to-noise map 
to have a pixel histogram with a standard deviation of 1. This 
process scales the noise map with a factor of 1.5. The reseated 
noise map is shown as contours in Fig. |2] together with the the 
signal-to-noise map which is constructed by dividing the signal 
map with the scaled noise map. 

The maps have a pixel size of 4". This means that the over- 
sampling factor compared with the original resolution of the ob- 
servations is 19.5/4 ~ 5. Such a fine pixelization is preferen- 
tial for the map-making process to be effective. The resulting 
map has a pixel-to-pixel noise level which is affected by high 
frequency noise on a scale smaller than the beam. This noise 
(which is due to small pixelization) can easily be removed by 
smoothing. We smoothed the final maps using a 10" Gaussian in 
order to remove that high-frequency noise component and pro- 
duce cleaner maps. 

3.2.2. Iterative mapping 

Some artefacts such as "sidelobes" around point sources are seen 
in the final map because of the filtering. In order to remove those 
artefacts, the entire Mini crush reduction was applied a second 
time, but instead of building a source model from the actual data 
being reduced, we used a source model based on the results of 
the first reduction. The part of the map with a signal above 4.5cr 
was used as the source model. 
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Fig. 3. Radial profiles of the brightest source (source #1): in 
the first iteration of the data reduction (solid line and boxes), 
a source model was successively built during the data reduction; 
in the second iteration (dashed line + circles), a 4.5o--clipped 
version of the first map was used as a source model. The itera- 
tive method makes it possible to recover flux that was lost at the 
first stage of the data reduction. 
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Fig. 4. Solid curve: Histogram of the pixel distribution in the 
central 10 arcminutes of the signal-to-noise map. Dashed curve: 
Histogram of the mean pixel distribution in the same region of 
the jackknife maps, as described in the text. The jackknife proce- 
dure has removed the excess positive signal (due to the sources) 
and the negative signal (due to the sidelobes) that are present in 
the signal map. 



Figure [3] illustrates the result of that iterative process: radial 
profiles of the brightest source in the map are shown for each 
of the two iterations. The difference between the two profiles is 
significant. We have observed that the difference is greater the 
brighter the source. The iterative process makes it possible to 
recover the flux that was lost in the first iteration. We have com- 
pared the fluxes of our two most significant sources (after two it- 
erations) with the values obtained using the standard data reduc- 
tion pipeline that was used for the calibrator, Uranus. The agree- 
ment showed that two iterations we re sufficient, as also found by 
other groups (e.g. lWeifi et al.ll2009l) . 

4. Results 

Figure |2] shows the signal-to-noise map obtained from the data 
reduction process described above. Several sources are visible. 
Since the noise level increases steeply toward the outer parts, 
sources with low signal-to-noise ratio in that part of the map 
may have high flux densities. In this section, we discuss the noise 
properties of the data and present the methods that we have used 
to identify sources and measure their properties. We also present 
Monte Carlo simulations performed to quantify the degree of 
completeness to which sources can be extracted, and to estimate 
the amount of flux boosting due to the confusion noise. Finally, 
we estimate the magnification of each source due to gravitational 
lensing, using a simple model of the Bullet Cluster and assuming 
that all sources are at a redshift of 2.5. 

As shown in the following analysis, neither completeness nor 
flux boosting corrections change the results very much, because 
of the conservative detection treshold that we have adopted. The 
results is a fairly robust catalog of 17 sources. 

4.1. Noise properties 

To quantify the noise level in the map, we constructed 500 so- 
called jackknife noise maps, which are obtained by multiply- 



ing half of scan maps (selected randomly) by (-1), before co- 
adding all the scan maps. The jackknife maps should there- 
fore be free from astronomical sources and from artefacts from 
the data reduction pipeline, and reveal the nature of the sta- 
tistical noise in the data. Each jackknife map was smoothed 
with a 10" Gaussian after coadding, giving a final resolution of 
V(19.5")2 + (10")2 ^ 22". Then, we calculated the mean of the 
500 histograms of the pixel distributions measured in the central 
10' of each jackknife map. 

Figure |4] shows the mean histogram, in units of signal-to- 
noise value, together with the histogram of the pixel distribution 
of the signal-to-noise map, extracted from the same central 10'. 
The jackknife histogram is fitted by a Gaussian with standard 
deviation 1 (as expected). If we instead make the same calcula- 
tion in the signal map, we find that the jackknife histogram is 
well fitted by a Gaussian function with mean /i = 1 .4 x lO"*" mJy 
and standard deviation cr - 1.17 mJy. Therefore, when exclud- 
ing astronomical sources and systematic effects from the data 
reduction, the statistical noise level is 1 .2 mJy/beam. 

4.2. Source extraction 

Although the noise in our final map is fairly uniform across the 
central 10 arcminutes, it increases slightly with radius across 
that area, and rapidly outside (see the noise map contours in 
Fig. 3). In order to identify significant sources in the map, a well- 
defined and mathematically j ustified algorithm must be used. 
Following other authors (e.g. iBeelen et al.ll2008h . we used the 
so -called "Gaussian ma tched filter" (GMF) technique outlined 
bv lSerieant et al.l(l2003h . This method is optimal for point source 
extraction in a sense, although the p erformance is degraded 
for crowded maps ( ISerieant et al.ll20()3l) . 

The GMF significance map, F/ AF, is computed as 

_ (S ■W)®P 
AF yfw~®P^ 
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Table 2. List of sources extracted from the LABOCA map. 



Source 


--.,/'Tonnn\ 
Q'(JZUUUJ 


Sil T'^AAA^ 


Flux density" 


Deboosted'^ 


Demagnified^ 


t j /\r 










flux density 


flux density 






(h:m:s) 


\ ■ ■ ) 


(mjy) 


(mJy) 


(mJy) 




1 


ri/'.CO.TT /lO 

Uo:5o:37.o2 


C C . C7 . A/f O 

-55:57:04.8 


48.6 ± 1.3'" 


48.0 +1.3 


0.64 


82.8 


2 


Uo:jo:24.4/ 


-jD:jj:12.3 


1 C 1 1 1 A 

13.1 ± l.U 


1/1-7 1 1 A 
14. / + l.U 


8.8 
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06:58:45.60 
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7.8+ 1.5 
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5.2 


11.9 


12 


06:58:52.22 


-55:55:45.7 


5.5+1.2 


4.5 + 1.2 


3.4 


11.2 


13 


06:58:22.88 


-56:00:40.7 


4.8+ 1.2 


3.8+1.3 


2.9 


11.0 


14 


06:58:46.68 


-56:02:11.8 


7.2+ 1.9 


4.6 + 2.5 


3.8 


10.8 


15 


06:58:33.69 


-55:54:40.8 


4.6+ 1.1 


3.6+1.2 
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10.1 
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06:58:12.44 


-55:57:29.7 
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1.9 


9.2 


17 


06:59:15.72 
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23.6 + 5.9 


f 
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Notes. Statistical uncertainties on the listed positions are 1-2", which is smaller than the pointing uncertainty. Flux density as extracted from 
the map. Flux density connected for boosting due to confusion noise. Flux density corrected for lensing. Significance of the detection in 
the Gaussian-matched-filtered map. Source #1 is extended relative to the 22" beam: it has an apparent size of 29.2" x 23.3". Source #17 lies 
in the outer part of the map where the noise level is high and the method used to deboost the flux densities fails. 



where S is the signal map, W is the weight map (the reciprocal 
of the noise map, squared), and P is a Gaussian of the same size 
as the beam. The i8> sign denotes a convolution. 

We generated two GMF maps, one for each iteration per- 
formed in Minicrush (see Sect. 3.2.2). Setting a threshold 
F/AF > 9, which roughly corresponds to a signal-to-noise level 
of 4, we extracted 19 sources in the first GMF map and 22 
sources in the second. We decided to be conservative and not 
include in the final catalog the sources that had appeared in the 
GMF map of the iterated map. In addition, we excluded the two 
sources that were present in the first GMF map but not in the 
second. Our final source list thus contains 17 sources. 

Figure |5] shows the GMF map calculated using the iterated 
signal map. The black contour indicates the 3 mJy/beam noise 
level, and the red circle has a diameter of 10 arcminutes. The 17 
identified sources are marked. 

Whereas source finding was done in the GMF map, measure- 
ment of the properties of the identified sources was done in the 
signal map. The following scheme was used: 

1 . We searched for the peak value in the GMF map and fitted a 
two-dimensional Gaussian. 

2. The position obtained from that fit was used as a starting 
guess in the fit to the real source in the signal map. The fit- 
ted function is the sum of a two-dimensional Gaussian and 
a tilted plane. The use of a plane reduces the effect of re- 
maining sidelobes around bright sources. Those are due to 
the filtering in the data reduction and are not completely 
removed by the iterative method described in Sect. 3.2.2 
(Fig E). Unless the source is clearly extended, the width of 
the Gaussian was fixed to that of the beam. 

3. The integrated flux density of each extracted source was cal- 
culated by integrating only the Gaussian part of the fitted 
function. 

4. The fitted Gaussian was then subtracted from the GMF map 
and a new search for the peak was performed. 



5. We iterated over the previous steps until the peak value in the 
GMF map was smaller than 9. This value was inferred from 
the simulations described in the next section. 

Table |2] lists the properties of the 17 sources, sorted by de- 
creasing value of F/ AF, given in the last column. The fitted posi- 
tions and integrated flux densities are also given. All sources but 
Source #1 are point sources. The uncertainties on the extracted 
flux densities are estimated from the 500 jackknife maps. 

4.3. Completeness 

In order to understand the systematics of the source extraction 
procedure and quantify the extent to which the values of AF/ F 
correspond to real sources, we turned to simulations. A Gaussian 
source of the size of the beam was added at a random location 
within the central lO'-diameter region of a randomly selected 
jackknife map, and the Gaussian filtered map was produced, us- 
ing Eq. ([T]). We stepped through a range of values of the flux den- 
sities of the simulated source, ranging from 1 to 15 mJy, with a 
spacing of 0.5 mJy. For each flux density value, 500 sources were 
simulated and placed at a random location, and the GMF maps 
were used to find the simulated sources. Once they were found, 
a two-dimensional circular Gaussian was fitted to the simulated 
signal map at the same location. For each flux density bin, we 
then extracted information about the completeness and the re- 
covered flux as a function of the input flux density value. 

Figure |6] shows the results. The flux boosting is calculated 
as the fraction of the measured flux of a source to the input flux 
density. At flux densities larger than 6 mJy we see that, in the ab- 
sence of confusion noise, no boosting is observed. For lower flux 
densities the flux boosting becomes larger than one, so that ex- 
tracted sources are those that are placed on positive noise peaks, 
making them rise above the noise. Sources randomly placed in 
noise voids are not extracted. 

The completeness is defined as the fraction of sources that 
are recovered from the simulated maps out of the 500 input 
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Fig. 5. Gaussian filtered map, with the 17 detected sources 
marked with circles. The numbering of the sources is the same 
as in Table |2] The black contour corresponds to the 2 mJy/beam 
level in the noise map and the dashed circle marks the central 
10' area of constant noise level used in the analysis. 
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sources. We see that at 6 mJy the observations are 100% com- 
plete. At the level of 4.6 mJy, which is the flux of the most dim 
source in our sample, the completeness is ~ 85%. Therefore, our 
sample is highly complete, and completeness corrections are not 
very important. 

Other groups use instead the final signal map when simulat- 
ing the completeness (eg. Scott et al. 2008), but we wanted to 
focus our attention in these simulation on the effect of statisti- 
cal noise in the map. No source confusion noise is present in 
the jackknife maps. In the following section, we discuss the ef- 
fect of confusion noise and use noise-free sky realizations with 
a Schechter distribution to quantify the flux boosting for an in- 
dividual source extracted from the real map. 



4.4. Flux boosting due to confusion noise 

Several SMG surveys have shown that th e number counts 
steepen towards h igher flux densities (e.g. IScott et aTl 120061 : 
ICoppin et al.ll2006h . Therefore there are many more sources at 
low than at high flux densities. Most of those faint sources are 
below the noise level of submm maps, but they influence pho- 
tometric measurements of extracted sources, acting as a "sea" 
of sources, often referred to as "confusion noi se". This effect 
has be en discusse d by e.g. [Con don ( 1974i) and iHogg & Turned 
(Il998l) . Recently, ICoppinet aP (l2005h discussed the eff'ect of 
confusion noise on flux boosting in the SCUBA Groth strip sur- 
vey, and used a Bayesian technique to "deboost" the fluxes. We 
employ a similar technique to estimate the amount of flux boost- 
ing for our detected sources, which we describe in Appendix lAl 
The derived flux boosting corrections are small for most sources. 
The deboosted flux densities are listed in Table |2] 



Fig. 6. Results of the Monte Carlo simulations performed to es- 
timate the degree of completeness of the source extraction algo- 
rithm and the accuracy of the measured flux density at a given 
value of the input flux density. 



4.5. Lensing correction 

We built a simple lensing model of the Bullet Cluster in order to 
estimate the magnification of our observed submm sources. The 
model consists of two spherically symmetric components repre- 
senting the main component of the Bullet Cluster and the sub- 
component (the actual "bullet" to the west). Figure|7](left) shows 
contours of the projected mass di stribution inferred from weak- 
lensing observations by IClowe et a l. (2 0061) . overlaid on our 
870 fim i mage. We used the fi ts to the weak-lensing observations 
made bv lBradac et al.l (l2006h to the masses within a certain pro- 
jected radius R of each component, M(< R) - Mo(R/250 kpc)", 
with Mo = 0.22 x IO'^Mq and n = 0.8 for the main cluster 
and Mo - 0.17 x IO'^Mq and n = 1.1 for the subcluster. To 
place the two mass comp onents, we us ed the information given 
in Table 2 of the paper bv lBradac et al.l The redshift of our simu- 
lated Bullet cluster was set to 0.296 and that of the source plane 
to z = 2.5. Note that the magnification values are not very sensi- 
tive to the redshift of the sources: varying it from z = 2 to z = 3 
changes the magnifications by less than 10%. The magnifica- 
tion map was calculated following the derivation in the book by 
Schneider et al. ( 1992). This numerical calculation provides re- 
sults along the lines of those ob t ained analytically for two point 
masses bv ISchneider & WeissI d 19861) and for two isothermal 
spheres bv lShin & EvansI (l2008l) . but in the case of two power- 
law projected mass distributions. Because the model does not 
include lensing by the individual cluster galaxies, the location of 
the critical lines dilTer from the observed ones by about 10 arc- 
sec. The true magnification for a given source must therefore 
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differ from our derived values. Nevertheless, this simple model 
makes it possible to estimate the individual magnifications and 
the average magnification in a certain region. Using the calcu- 
lated magnifications, we corrected the measured flux densities. 
The lensing-corrected values are listed in Table |2] 



5. Discussion 

5.1. Searching for infrared counterparts to tine submm 
sources 

The identification of optical counterparts to SMGs is often diffi- 
cult due to their optical faintness, as much of the starburst lumi- 
nosity is highly obscured by dusty clouds. 

Rest-frame infrared photometry of SMGs is usually possi- 
ble because the extinction is much smaller than in the optical. 
With the Sp itzer satellite an d its imaging pho tometers IRAC 
jFazio et al.ll2004 and MIPS jRieke et al.ll2004l) . it is now pos- 
sible to obtain high-resolution infrared images of these distant 
objects. The number counts of sources at 3.6 /im is high and 
counterpart identification is complicated by the large positional 
uncertainty of the submm source (see below). In the 24 jum band, 
the number density of sources is smaller and a more secure coun- 
terpart identification can be made, if a source can be identified. 

Figure[8]shows cutouts of IRAC 4.5 and IRAC 8.0 and MIPS 
24.0 fim images of 16 sources detected in the LABOCA map. 
Source #17 has only partial Spitzer coverage and is omitted from 
this figure. The Spitzer data used in this study are described in 
Sect. 12.21 and the coverage with respect to the LABOCA map is 
shown in Fig.|7] 

The full I'xl'cutouts are only shown for reference here. It is 
not necessary to search for infrared counterparts in such a large 
field. The region of interest (within a certain search radius) is 
the sum of the contributions from the mean pointing error of 
LABOCA of ~ 4"; the statistical uncertainties on the fitted po- 
sitions, which are 1 - 2"; the systematic uncertainty due to con- 
fusion noise, ~ 3 - 4" and a possible misalignment between the 
LABOCA and Spitzer maps of the order of 5". Together, this 
adds up to a search radius of ~ 10". 

Table [3] lists the coordinates of li kely Spitzer counterparts , 
found using the software SExtractor (iBertin & Arnoutslll996h . 
and their measured flux densities in the 3.6, 4.5, 5.8, 8.0, and 
24 i^m bands. We extracted sources with six adjacent pixels 
above 3.5cr, and used an aperture of five pixels for the photom- 
etry; applying the aperture corrections listed in the IRAC and 
MIPS data handbooks. The listed uncertainties on the flux den- 
sities are the statistical errors given by SExtractor We estimate 
the systematic uncertainties to be ~ 10%. In the cases where a 
source has been extracted in one IRAC channel and not another 
we list upper limits. 

We found in total 9 sources with infrared counterparts, and 
two where counterpart identification is complicated by the larger 
number of sources in the short wavelenght IRAC bands. For 
two sources we see an excess of flux in the MIPS map, but it 
is not significant enough to be extracted. Deeper MIPS imag- 
ing would therefore be useful. Deep high-resolution radio maps 
could also be used to distinguish between sources in the shorter 
wavelength bands. Comments about the individual sources are 
presented below. The lack of infrared counterparts was to be 
expect ed: for comparison, in the SHADES survey iPope et all 
(I2006h found 21 secure Spitzer counterparts to the 35 submm 
sources. An important result is that the large positional uncer- 
tainty of LABOCA (~ 10" for most sources) is reduced because 
the pointing accuracy of Spitzer is less than 1". 



Out of the nine sources with likely infrared counterparts, five 
have sufficient coverage to make it possible to investigate the 
shape of the mid-infrared SED. Ivison et al. (2004) suggested a 
diagram based on mid-infrared colors 5'8.o/5'4.5 versus 524/5' g.o- 
Based on the redshift tracks (the position as a function of red- 
shift for a SED in color-color space) of typical starburst and 
AGN-type spectral energy distributions, such a diagram could 
distinguish between strong starburst S EDs and powerful A GN. 
A similar diagram was a lso used by Ivison et al.l ( l2007h and 
iBeelen et aP (l2008l) . while lHainUne et al.l(l2009l) showed that its 
diagnostic capacity is limited. The five sources with identified 
counterparts in these three Spitzer bands have colors that lie in 
the starburst part of the diagram. This does not exclude the pos- 
sibility of contributions from AGN to the dust heating, but indi- 
cates that the galaxies are starburst, rather than AGN, dominated. 
Infrared spectroscopic measurements could be used to investi- 
gate further the power source of those galaxies. 



5.2. Notes on individual sources 

Aside from Source #1, none of the other sources have been 
detected previously in the mm or submm. Our observation of 
Source #1 is discussed in the context of other observations. The 
few other sources for which complementary observations exist 
are discussed as well. 

Source #1: With a deboosted flux density of 48.0 + 1 .3 mJy, 
this source is one of the brightest SMGs ever detected around 
870 ijm. This is very likely because of its proximity to a caus- 
tic line , which provides a larg e magnification. From their lensing 
model. lOonzalez et al.l (l2009l) estimated that the two brighter im- 
ages of the galaxy, A and B, which are separated by 8.6", have a 
magnification of 25 and 50, whereas the third image, C, located 
47" away from the first one, would have a magnification four 
times lower than the first one. Our final map was smoothed to 
22" and our fit to the Source #1 gives a position between the im- 
ages A and B; the flux that we measure comes most likely from 
those two images. The total magnification is therefore 75. 

In the LABOCA study of the protocluster J2 142-4423 by 
iBeelen et al.l (l2008l) the brightest source has a flux density of 

21. 1 mJy. The flux density at 850 fim, which is one of the 
SCUBA wavelengths, can be extrapolated using a submm spec- 
tral index of 2.7, giving a flux density 7% higher than at 870 //m. 
The SCUBA flux of Source #1 would therefore be around 
51 mJy. For c omparison, the brightest source in the SHADES 
survey (ICoppin et al. 2006) has 22 mJy (and a signal-to-noise 
rati o of 4.9). In the submm survey of massive clusters of galax- 
ies, iKnudsen et al.l (l2008h detected two bright sources towards 
Abell 478 and Abell 2204, with flux densities of 25.0 and 

22.2 mJy respectively. 

Source #1 was detected using the AzTEC bolometer array at 
1.1 mm (W08); we can compare the astrometric position, angu- 
lar size and integrated flux density with the values measured at 
870 fim. 

- Position: the distance between the estimated central position 
of the two sources is 4.1". This is well within the pointing 
error margins of the two telescopes. 

- Source size: W08 reports a source size of 36 ± 1.3" x 
32 ± 1.2". To compare that with the LABOCA source, we 
smoothed our map to the AzTEC resolution of 30", and fit- 
ted a two-dimensional Gaussian. The fitted FWHMs at this 
resolution are 36.2 + 1.3" x 30.1 + 1.2", in good agreement 
with the value of W08. 
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Right Ascension (J2000) Right Ascension (J2000) 



Fig. 7. Left: Signal-to-noise map with overlaid circles indicating extracted su b mm sources. The contours show the 
the projected mass density from the weak lensing analysis by iClowe et al.l (l2006l) . The contours range from 40 
to 85% of the maximum value and are spaced by 15%. The weak lensing map was retrieved from the website 
http : //flamingos . astro .ufl . edu/le0657/public .html. The rectangles show the regions of complete coverage of Spitzer 
MIPS (small rectangle) and IRAC (large rectangle). Right: Signal map (in units of Jy/beam) with contours of the X-ray surface 
brightness from XMM-Newton observations. The noise level in the signal map increases rapidly towards the outskirts because of 
the low coverage there. 

Table 3. Photometry of possible infrared counterparts to the LABOCA sources. 



Source 


a(y2000) 


5(72000) 


^laboca" 


S 3.6^m 


S 4.5 




S 8.0 ;im 






(h:m:s) 


r :':") 


(") 


(A'Jy) 


(A-Jy) 


(A-Jy) 


O^Jy) 


(A'Jy) 


2 


06:58:24.59 


-55:55:12.2 


1.1 


71.2 ±0.3 


99.5 ±0.5 


77.8 ± 1.6 


53.1 ± 1.9 


507.5 ± 20.5 


5 


06:58:27.17 


-56:01:16.8 


1.0 


12.1 ±0.2 


15.1 ±0.3 


< 19.9 


24.9 ± 1.5 


249.6 ± 24.3 


6 


06:58:28.92 


-55:53:52.0 


3.6 


13.7 ±0.4 


8.8 ±0.5 


<6.0 


< 2.0 


<40 


7 


06:59:00.72 


-55:52:21.9 


6.8 


89.2 ±0.2 


54.6 ± 0.4 


26.8 ± 1.5 


< 18.7 


<40 


8 


06:58:23.96 


-55:57:19.5 


3.6 


55.7 ±0.6 


68.5 ±0.6 


50.1 ± 1.5 


40.5 ± 1.5 


582.3 ± 24.4 


10 


06:58:45.39 


-55:58:46.4 


2.3 


17.5 ±0.7 


23.7 ±0.7 


28.2 ± 1.9 


< 18.4 


608.0 ± 39.0 


11 


06:58:53.52 


-56:00:48.4 


4.3 


33.4 ±0.4 


47.0 ± 0.4 


72.2 ± 2.2 


88.5 ± 1.6 


228.8 ±33.5 


15 


06:58:34.14 


-55:54:37.2 


5.3 


46.1 ±0.3 


37.2 ±0.4 


37.4 ± 1.4 


30.1 ±2.1 


301.3 ± 15.5 



Notes. Distance between the IRAC 1 position and the central LABOCA position. 



- Flux density: After removal of the contribution from the 
Sunyaev-Zeldovich effect, the flux density of the sou rce at 
1.1 mm was estimated to 13.5 mJy (IWilson et al.ll2008h . This 
value, however, is too low and is being revised to about 
20 mJy (G. Wilson, priv. comm.), which would give a spec- 
tral index a^^y^ -3.7 between the two measurements. 

Far-infrared observations of this source by the BL AST ex- 
periment were recently published bv lRex et al.l (l2009h . BLAST 
observed at wavelengths 250, 350 and 500 fim, with FWHM 
beam sizes of 36, 42 and 60", respectively. High-significance 
detections were reported in all three bands; the flux measure- 
ments were contaminated by that of an elliptical galaxy, which 
is a spectroscopically confirmed member of the Bullet Cluster 
and which had to be modeled and subtracted. After making cor- 
rection for the SZ flux and color (due to the width of the BLAST 
filters), the flux densities in the three bands are 94 + 30, 96 + 27 
and 110±21 mJy for bands ranging from shorter to longer wave- 
lengths. 

Figure |9] shows the spectral energy distribution of 
Source #1. The measured data points are from from Spitzer, 



BLAST, LABOCA a nd AzTEC. We compared the models of 
iLagache et al.l (|2003|) of template starburst galaxies to those 
data points. We used a magnification of 75, assuming t hat the 
flux co mes from images A and B discussed by Gonzal ez et al.l 
(120091) . We display the SED of a starburst g a laxy a t redshift 2.9, 
which is the redshift estimate of iRex et al.l (l2009l) . with a total 
luminosity of 10" '^L©. The SED of a more redshifted galaxy, 
at redshift 3.9, with a luminosity of 1O"'^L0 is in better agree- 
ment with the long-wavelength data. However, it does not fit the 
24 fim value. Although we have not done a formal fit, it seems 
that the current data would require a different template SED. 
Alternatively, the uncertainties on the data points may have been 
seriously underestimated. 

We do not list any Spitzer photometry of source #1 because 
accurate measurement would require not only processing the raw 
data but also a careful subtraction of the foreground elliptical 
that lies between the t wo images of that sou rce. Such a work 
has been performed bv lGonzalez et al.l ( l2009l) and we used their 
quoted values in the analysis. 
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Fig. 9. Spectral energy distribution of Source #1. The curves 
show modeled SEDs of starburst galaxies taken from the tem- 
plates of Lagache et al. (2003) magnified by a factor of 75. The 
magnification factor is taken from the lensing model used by 
Gonzalez et al. (2009) who find a magnification of 25 for the 
first image and of 50 for the second. The dotted-dashed curve 
corresponds to a starburst galaxy of intrinsic total luminosity 
Ltot - 1O"'^L0 redshifted to z = 2.9. The solid curve shows 
the SED of a starburst galaxy of Ltot = lO" ''L© redshifted to 
z = 3.9. The data poin ts are the 24 //m Spitzer measurement of 
iGonzalez et aT ('2 009h. t he 250 , 350 and 500 yum BLAST mea- 
surements of Rex et al.l (l2009h . our 870 yL/m LABOCA point, 
and the AzTEC 1 . 1 mm measurement of W08. 



Source #2: This is our second most significant detection at 
Gaussian filter value 29.9 or signal-to-noise ratio S/N^15. No 
previous mm or submm detection of this source has been re- 
ported, although in the map showing the detection of Source #1 
in W08 (their Fig. 1) there is an indication that they also detected 
this source. 

Fig. [8] shows that this galaxy has likely counterparts in all 
available Spitzer bands. There is a source to the north of the 
center of the LABOCA detection in the 3.6 and 4.5 yum maps 
that might be the reason for the apparent elongation of the red 
contour towards the north. However, this source is outside the 
10" search radius and does not show up in the 24 yLVS\ map. 

Source #3: A bright, extended Spitzer source is seen to the 
southeast of the LABOCA position. Its large angular size at 
8.0 yum and its distance from the submm position indicates that 
it is not the infrared counterpart. It is detected in the 2MASS 
catalog, but no redshift is indicated. Two sources are detected 
within the 10" circle. The source closest to the LABOCA posi- 
tion is identified in all IRAC band, whereas the other source is 
only found in IRACl and IRAC2. Both these sources are very 
faint due to this source confusion we choose to not list a coun- 
terpart. A flux excess is identified in the MIPS map, but it is not 
significant enough to be extracted by SExtractor. 

Source #4: The identification of a counterpart to the submm 
source is complicated by the bright foreground source seen in all 
five Spitzer bands. This source, at z = 0.097, is also found in 
the 2MASS catalog. In radio observations with the ATCA array, 
iLiang et al] (|200C^ detected a source at the position of the fore- 
ground galaxy. The flux density at 1.344 GHz is 0.8 ± 0.05 mJy 



and drops towards higher frequencies. The source is not detected 
at 4.4 and 8.8 GHz. 

Due to the positional offset, it is not very likely that the 
submm emission seen in the LABOCA map is caused by the 
z - 0.097 galaxy south of the submm detection. There is a possi- 
bility that a small part of the total LABOCA flux density is com- 
ing from this foreground source, but not a substantial fraction. A 
faint source is detected to the north of the central position, but 
its position is slightly outside the 10" search radius. 

Source #5: Two sources with detections in the IRAC 1, 2 
and 4 bands are found close to the central position. A MIPS de- 
tection is associated with the centralmost source, which is situ- 
ated only 1" from the LABOCA position. The other source lies 
6" away which point towards a counterpart association with the 
central source. High resolution radio imaging is needed to com- 
pletely differentiate between these two sources, but we list the 
most central source here as a counterpart. Another source, to the 
south-east, is detected only in IRACl, and is likely too faint to 
be a counterpart. 

Source #6: A source is detected in the IRAC 1 and 2 bands 
south of the LABOCA position. The brighter source north-west 
of the central position is further out than the search radius. 

Source #7: A significant source is detected in the IRAC 1, 2 
and 3 bands west of the center of this source. Its SED is declining 
towards the longer wavelength IRAC bands. 

Source #8: A counterpart is identified in all Spitzer bands to 
the north of the LABOCA central position. 

Source #9: Almost 10" from the central position lies a sig- 
nificant source, which is probably too far out and too bright to be 
the SMG counterpart. Closer to the central position two IRACl 
sources are extracted, but their similar flux levels and distance 
to the center makes the counterpart association ambiguous. We 
therefore choose to not list a counterpart for this SMG. 

Source #10: A bright MIPS source with counterparts in 
IRAC 1, 2 and 3 is situated close to the LABOCA central po- 
sition. This is a tentative counterpart to the LABOCA source. 

Source #11: A possible counterpart is detected in all the 
Spitzer bands. 

Source #15: The source to the northeast is identified as a 
counterpart and exists in all the Spitzer bands. 

Source #17: lies outside the R - 5' circle from the center 
in a region when the noise is high. Our deboosting algorithm 
failed in that case, so we only quote the measured flux density 
of 23.6 + 5.9 mJy. This very bright source is also present in the 
AzTEC 1.1 mm map (D. Hughes & 1. Aretxaga, priv. comm.). 

5.3. Number counts 

Figure [TO] shows the cumulative number counts, defined as the 
number density of sources with flux density larger than S , de- 
noted A^(> S), derived from our observations and from other 
studies. There is good overall agreement between our results and 
those from previous LABOCA and SCUBA surveys, although 
our surveyed area is not as large as that of other studies, as re- 
flected by the larger errorbars. 

Because of the gravitational magnification, the surveyed area 
in the source plane is smaller than that in the image plane. We 
used the simple lensing model described in Sect. 4.5 to calculate 
the mean magnification in the central part of the map used for the 
analysis. We calculated a mean magnification factor of 3.3, i.e. 
the area of the source plane is 3.3 times smaller than the surveyed 
area in the image plane. The magnification of the source fluxes 
is also taken into account by using the de-magnified source flux 
densities listed in Table |2l 
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Table 4. Cumulative number counts based on the central 10 ar- 
cminutes of our map. 
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Fig. 10. Cumulative number counts derived from the central 10 
arcminutes of the 870 fim map, after lensing and flux boost- 
ing correction. As the observations are highly complete, no 
completeness correction has been made (see the discussion in 
Sect. 14.31 and Fig. |6]l. For comparison , the result i ng nu mber 
counts from the LA BOCA surveys of IWeifi et all (|2009|) and 
iBeelen et all (l2008h. and the SCU BA survevs of ICoppin et alJ 
(l2006l) and iKnudsen et al.l (l2008h are displayed. The arrow on 
our point in the lowest bin indicates a lower limit: we reach the 
sub-mJy regime only because of the one bright source with an 
estimated magnification of 75. The surveyed are for this mag- 
nification is much smaller than the mean magnification factor 
of 3.3 described in the text, therefore this point only denotes a 
lower limit. 



Our counts were inferred from the 1 3 sources found within 
the central 10 arcminutes (the region within the red circle shown 
in Fig. Is]), where the noise is uniform. From our 13 sources we 
construct the number counts in four flux-bins. We probe the 
counts down to ~ 0.64 mJy, a regime which can be investi- 
gated only by utilizing foreground lensing clusters. The number 
of points per flux-bin and the counts are displayed in TableU] We 
estimate errorbars for each data point from Poi ssonian statistics, 
using the tabulated values from lGehrelsl (Il986l) . 

In Fig. [To] we also show the results from the LESS survey 
(filled diamonds), which covered an area of 30' x 30' of the 
Extended Chandra Deep Field South to a uniform noise level 
of 1.2 mJy/beam dWeiB et al.l 12009.) . That survey has a noise 
level comparable to ours, but it covers an area which is 10 times 
larger The apparent deficit of galaxies in t he LESS has also 
been demonstrated in other wavebands (see IWeifi et aP (l2009t) 
and references therein). The results of the LABOCA stud y of 
the protocluster J2 142-4423 at z = 2.38 dBeelen et alJl2008h are 
plotted in Fig.[lO]as open circles. The excess seen at 5 > 3 mJy 
might be an effect of clustering of submm sources in the proto- 
cluster. 

We also c ompa re with the SCUBA surveys from 
(120061) (SHADES, filled circles) and 
(I2008h (open boxes). The SHADES survey. 



Copp in et alj_ 
Knudsen et al 



which covered blank fields, probes the high-end of the number 
counts, while the Kniidsen et alj survey targeted 10 lensing 
galaxy clusters and reached lower flux levels. A lt hough our 
surveyed area is similar to that of IKnudsen et all (l2008l) . we 
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find only one sub-mJy source compared to the seven found by 
IKnudsen et al.l This is because their survey covers a larger area 
of high magnification. 

5.4. Resolving the cosmic infrared background 

It is interesting to investigate how much of the cosmic infrared 
background (CIB) flux is resolved in our LABOCA observation. 
The CIB was discovered by the COBE iri struments FIRAS and 
DIRBE (iFixsen et al.lll998l : iHauser et al.l [1998). Those experi- 
ments detected an isotropic infrared background signal that was 
thought to originate from the integrated effect o f star formation 
in the history of the universe (iDwek etal.lll998h . But the coarse 
resolution of the COBE satellite made it impossible to resolve 
the sources responsible for the background signal. The detection 
of those sources had to wait for the invention of mm and submm 
bolometer cameras on large ground-based telescopes, resulting 
in angular resolutions of tens of arcseconds. 

Because of the steep number counts, most of the ffux 
coming from submm galaxies originates in low-flux sources. 
Therefore, observations of lensing foreground clusters that probe 
the faintest number counts potentially will resolve a large frac- 
tion of the CIB. 

We use the central region of 10' diameter for this calculation. 
The amo unt of CIB flux in th at area, at 870 yum is estimated from 
Eq. I of lDweket all (Il998h and yields ~ 900 mJy. The total 
flux in the sources extracted from the LABOCA map within that 
region (the same sources that are used for the number counts, but 
without the flux correction due to magnification) is 124+4.0 mJy. 
This calculation is valid because gravitational lensing preserves 
surface brightness. Thus, 14% of the CIB has been resolved in 
the LABOCA observations. 



5.5. Residual SZ emission 

Although the data have been filtered to remove the extended 
emission, some level of residual SZ emission may be present. 
Close inspection of the central region of the map shows extended 
emission between the two brighter sources. Source #1 and #2. 
We note that W08 (their Fig. 1) see a similar extended structure 
around the same central sources. In order to set an upper limit 
to the amount of residual SZ signal present in the region of the 
map covered by the Bullet cluster (which corresponds roughly 
to the area of a disk with an outer radius of about 2 arcminutes), 
we calculated the difference between the excess flux density in 
that region (estimated from Fig. [U and the flux density in re- 
solved point sources in the same region. Removing the positive 
noise contribution estimated from the jackknife noise maps we 
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find that 40 mJy residual SZ emission might be present in the 
final map. 

We can use the result s of the APEX-SZ obse rvations of the 
decrement at 150 GHz bv lHalverson et al] (l2009l) to estimate the 
SZ flux density across the sky area covered by the Bullet Cluster. 
iHalverson et al.l (1200^ fitted an elliptical beta model using an X- 
ray prior on j8 = 1.04^°' [[] and found a core radius 0^. = 142 + 18", 
an axial ratio of 0.89 ± 0.07 and a central Compton parameter 
yo - 3.31 + 0.30 x 10"^. Integrating the beta model over the 
central region with a maximum radius of 2 arcminutes, and tak- 
ing for simplicity a spherical model with the corresponding core 
radius, we estimate that the total SZ flux density from that area 
is of the order of 400 mJy. Therefore, we cannot exclude that a 
maximum of 10% of the total SZ flux of the Bullet Cluster may 
remain in the map in the form of extended emission. 

6. Conclusions 

Continuum observations at 870 yum of the Bullet Cluster have 
been presented. The data were filtered to remove large-scale sig- 
nal such as the Sunyaev-Zeldovich increment from the cluster. 
The main results are summarized below: 

- Seventeen submm sources with signal-to-noise ratios larger 
than 4 were detected in the map. Their measured fluxes den- 
sities range from 4.6 to 48 mJy. For each source but one, 
which lies in a noisy area, we calculated the value of the flux 
density corrected for the flux boosting due to confusion noise 
and for the gravitational magnification by the cluster 

- The brightest submm source coincides with a previously 
reported galaxy at an estimated redshift of 2.7-2.9 de- 
tected by Spitzer and the AzTEC and BLAST experiments. 
With its flux density of 48 mJy it is one of the brightest 
submm galaxy ever detected. After correction for gravita- 
tional magnification |yu|, the intrinsic flux of the source is 
about 0.64 (H/75r'mJy. 

- We found reliable infrared counterparts for nine of the 
submm sources. An infrared color-color analysis suggests 
that they have starburst-dominated spectral energy distribu- 
tions. 

- The cumulative number counts derived from the observa- 
tions agree well with those from other surveys. 

- The observations resolve 14% of the cosmic infrared back- 
ground radiation at 870 /im, in a sky area of ~ 78.5 arcmin^. 

We plan to apply a similar analysis to several other cluster 
fields for which we have LABOCA data. The results presented 
here will also be used to remove the contribution of the submm 
sources and obtain a map of the SZ increment at 870 /im in the 
Bullet Cluster with sub-arcminute angular resolution. Although 
the SZ decrement of the Bullet Cl uster has been observ ed using 
several instruments (e.g SES T bv lAndreani et al.ll 19991 AP EX- 
SZ by IHalverson et all l2009l ACT by 'Hincks et al."2009', and 
recently the South Pole Telescope by Plagg e et al. 2009) the 
only observation of the increment was do ne using ACBAR a t 
275 GHz, with an angular resolution of 4.5' dGomez et alJ2004l) . 
Removal of lensed background sour ces is not a trivial task, si nce 
it can also bias the SZ measurement dLoeb & Refregiedl997h . In 
the case of the Bullet Cluster, the extremely bright lensed submm 
source close to the center of the cluster field has an integrated 
flux density comparable to that of the SZ flux from the central 
region of the cluster. A careful analysis must be performed to 
recover the SZ increment in such a system, and may require a 
joint analysis using measurements in other wavebands. 
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Appendix A: Flux deboosting 

This appendix describes the flux deboo sting algori t hm us ed in 
this work, following the prescription in IScott et aP (l2008l Sec. 
5.1), whose notation we use in the rest of this section. 

The analysis builds upon Bayes theorem. It gives, for a 
source with a measured flux density S,,, and uncertainty cr,„, the 
posterior probability distribution of its intrinsic flux density S i: 

P{S i \S,n, cr,„) = — , (A. 1 ) 

P{S,„,a,n) 

where P(S ,) is the probability distribution of intrinsic flux den- 
sities, P{S,„, <T,„ \S i) is the likelihood of the observed data, and 
P{S,„, cr,„) is a normalization. 

The likelihood of observing the data, P{S,„, cr,„ \S i), is mod- 
eled by a Gaussian distribution. This is valid because of the 
Gaussian distribution of pixel values in the jackknife maps (see 
Fig. ID . The UkeUhood is 

P{S„,(T,„ I 5,) = exp p^^'f'^' ) (A.2) 

V2^ \ 2o-2 / 

We calculated the prior sky distribution P(S i) from Monte Carlo 
simulations. We used a iSchechter Cl976i) form of the number 
counts 

dN I S \"*^ 

-=«'(-) exp(-5/Sa (A.3) 

as fitte d to the soti r ce cou nts in the SCUBA SHADES sur- 
vey of ICoppin etaP (l2006l) . We used N' = 1703deg Vjy 
S' = 3.1 mJy and a = -2.0, which have been scaled from the 
ICoppin et al.l values using a spectral index of 2.7. We then drew 
randomly source fluxes from that distribution and placed point 
sources (convolved with the beam) on noiseless sky maps. The 
position of each source was drawn from a uniform distribution, 
so no clustering was introduced in the simulation. We simulated 
a region of 10' x 10'. In order to reduce edge effects, each simu- 
lated map was 3' larger, and a smaller map was extracted. 

We generated ~ 10^ sky maps and measured the pixel distri- 
bution in each of them. The mean distribution of pixel values in 
these maps is a measure of P{S i). 

From those simulated sky maps, we could calculate the pos- 
terior flux distribution of a source extracted from the map with 
flux density S,„ and uncertainty cr„,: it is obtained by multiplying 
the prior P{S ,) (from the simulations) with the Gaussian distri- 
bution (eq. lA.2l ). and dividing by P{S,„, o"m)- 

Figure lATI shows examples of posterior flux distributions for 
three combinations of S,„ and cr„,. The amount of flux boosting 
for a certain source depends both on the signal-to-noise ratio and 
the value of the flux density; for example, for a signal-to-noise 
ratio of 4, a source with a measured flux density of 10 mJy will 
have been boosted more than a source with a flux density of 
5 mJy. 
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Fig. 8. I'xl'cutouts of the Spitzer 4.5, 8.0 and 24 fim maps centered on the positions of the LABOCA sources. The contours 
correspond to the 3.5, 4.5, 5.5 and 6.5cr-levels in the LABOCA map. The grey-scales for the Spitzer images are spaced between 
[0.0,0.4], [2.4,2.8] and [20.25,20.65] MJy/sr respectively and the pixel sizes are 0.85, 0.85 and 3.5 arcseconds. We have also 
marked the positions of the likely infrared counterparts to the submm sources. The counterpart association is described in detail in 
Sect.lOandlO 
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Fig.A.l. Solid curves: Posterior flux distributions for measured flux densities of 3, 5 and 7 mJy and a noise level of 1.2 mJy. 
Dashed curves: Measured flux distributions, modelled as Gaussian functions. Dotted-dashed curve: Prior flux density distributions 
estimated from noise-free simulations, populating maps with a Schechter flux distribution (Eq. ( IA.3I )). The prior (red dotted-dashed) 
is identical in all three panels. The posterior flux distribution is the result of multiplying the blue dashed curve by the red dot- 
dashed curve and normalizing. The two vertical lines in each panel show the measured flux (rightmost line) and the maximum of 
the posterior flux distribution (leftmost line). The difference between the two lines indicates the amount of flux boosting. 
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